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PREFACE 


This  report  is  part  of  RAND's  continuing  interest  in  the 
detection  of  high-altitude  nuclear  explosions.  In  particular,  it 
is  concerned  with  low-frequency  hydromagnetic  waves  similar  to  those 
generated  by  such  explosions,  and  which  are  detected  by  ground- 
current  measurements  at  world-wide  locations.  In  this  report,  we 
direct  our  attention  to  the  response  of  the  ionosphere  to  such  low- 
frequency  hydromagnetic  disturbances.  The  purpose  is  to  separate 
the  effects  of  the  medium  from  those  of  the  source  on  the  trans¬ 
mitted  disturbance. 
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SUMMARY 


In  this  paper  we  consider  the  propagation  of  hydromagnetic 
waves  at  altitudes  below  about  2000  kin.  This  work  is  an  extension 
of  a  previous  paper  (Greifinger  and  Greifinger,  1965)  where  heights 
below  about  500  km  were  considered.  We  treat  the  case  of  plane 
wave  propagation  in  the  vertical  direction,  and  assume  the  geomagnetic 
field  also  to  be  vertical  (polar  propagation) .  Currently  used  models 
indicate  that  the  ionosphere  (the  height  region  between  about  80  and 
500  km)  may  be  reasonably  represented  by  a  constant  Alfven  speed 
and  locally  exponential  ion-neutral  collision  frequency,  while  the 
lower  exosphere  (the  region  between  about  500  and  2000  km)  can  be 
adequately  described  by  an  Alfven  speed  which  increases  exponentially 
with  height.  By  using  these  approximate  forms,  the  solutions  of  the 
relevant  forms  of  the  electromagnetic  wave  equation  are  expressed  in 
terms  of  known  functions.  Analytic  expressions  for  the  magnetic 
transmission  and  reflection  coefficients  are  derived  and  analyzed, 
and  numerical  results  are  obtained.  These  results,  although  strictly 
applicable  only  to  high  magnetic  latitudes,  compare  favorably  with 
many  geographically  widespread  experimental  data.  Of  particular 
interest  is  the  prediction  of  a  prominent  double  transmission  resonance 
in  the  daytime  and  a  single  strong  resonance  at  night.  This  agrees 
with  the  measurements  of  Maple  (1959)  and  others.  Many  lesser  resonances 
are  also  found.  The  advantage  of  the  analytic  representation  is  the 
ease  of  interpretation  of  the  physical  results.  For  example,  simple 
expressions  describing  the  transmission  resonances  found  by  Jacobs 
and  Watanabe  (1962)  for  a  grounded  exosphere  are  found  to  be  a 


limiting  case  of  the  equations  derived  here. 
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FIGURE  CAPTIOi,  > 


Fig.  1 
Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Schematic  model  of  the  ionosphere  and  lower  exosphere. 

Magnitude  of  the  composite  magnetic  transmission  co¬ 
efficient,  |Tgl,  for  the  L-mode,  as  a  function  of  frequency 
for  daytime  sunspot  maximum  conditions .  For  purposes  of 
comparison,  the  transmission  coefficient  for  the  grounded 
exosphere  model  of  Jacobs  and  Watanabe.  and  the  ionospheric 
transmission  coefficient,  | tg |  ,  are  also  shown. 

Magnitude  of  the  composite  magnetic  transmission  coefficient, 
|Tjj|,  and  of  the  ionospheric  transmission  coefficient,  |tg|, 
for  the  L-mode,  as  a  function  of  frequency  for  daytime  sun¬ 
spot*  minimum  conditions.  In  this  figure  and  in  Fig.  2, 
which  both  correspond  to  daytime  conditions,  the  composite 
trar  nission  coefficient  shows  a  prominent  double  resonance. 

Magnitude  of  the  composite  magnetic  transmission  coefficient, 
|Tb| ,  for  the  L-mode,  as  a  function  of  frequency  for  night¬ 
time  sunspot  maximum  conditions.  For  purposes  of  com¬ 
parison,  the  composite  transmission  coefficient  for  a  uni¬ 
form  exosphere,  Tg(X  M  ») | ,  and  the  ionospheric  transmission 
coefficient,  |tg  ,  are  also  shown. 

Magnitude  of  the  composite  magnetic  transmission  coefficient, 
|Tg|,  for  the  L-mode,  as  a  function  of  frequency  for  night¬ 
time  sunspot  minimum  conditions .  For  purposes  of  com¬ 
parison,  the  composite  transmission  coefficient  for  a  uniform 
exosphere,  |T«(\  *  ») j ,  and  the  Ionospheric  transmission 
coefficient,  |tg|,  are  also  shown.  This  figure  and  Fig.  4 
show  that  the  ionospheric  resonance  Is  virtually  eliminated 
under  nighttime  conditions,  with  only  the  exospheric 
resonances  remaining. 

Magnitude  of  the  composite  magnetic  transmission  coefficient, 
|Tg|,  and  of  the  ionospheric  transmission  coefficient,  |  tg  | , 
for  the  R-mode,  as  a  function  of  frequency  for  daytime  sun¬ 
spot  minimum  conditions .  The  magnitudes  of  the  resonant 
peaks  are  much  smaller  than  those  shown  in  Fig.  3  for  the 
L-mode . 


Fig.  7  -  Magnitude  of  the  composite  magnetic  transmission  coefficient, 
|Tgj,  and  of  the  ionospheric  transmission  coefficient,  |tg|, 
for  the  R-mode,  as  a  function  of  frequency  for  nighttime 
sunspot  minimum  conditions .  The  magnitudes  of  the  resonant 
peaks  are  comparable,  in  this  case,  to  those  shown  in  Fig. 

5  for  the  L-mode. 


Fig.  8  -  Magnitude  of  the  composite  magnetic  transmission  coefficient 
for  the  grounded  ionosphere,  )Tg(H  *  0) | ,  for  the  L-mode, 
as  a  function  of  frequency  for  daytime  sunspot  maximum  condi¬ 
tions.  Grounding  the  ionosphere  has  caused  the  lowest  day¬ 
time  resonance,  shown  in  Fig.  2,  to  be  depressed  and  shifted 
upward  in  frequency. 
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Fig.  9  -  Magnitude  of  the  composite  reflection  coefficient,  |R|, 

and  of  the  ionospheric  reflection  coefficient,  |r|,  for  the 
L-mode,  as  a  function  of  frequency  for  daytime  sunspot 
maximum  conditions. 
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I.  INTRODUCTION 

One  o£  the  most  interesting  characteristics  of  both  natural 
and  artificially  produced  geomagnetic  micropulsations  is  the  existence 
of  narrow  bands  of  oscillations  (e.g. ,  Maple,  1959,  Tepley,  1961),  or 
equivalently,  the  existence  of  resonant  peaks  in  their  observed 
energy  spectra  (Ness  et  al. ,  1962;  Davidson.  1964,  Santirocco  and 
Parker,  1963;  Smith,  Provazeck  and  Bostick,  1961).  The  micropulsation 
energy  spectra  depend,  of  course,  upon  the  nature  of  the  sources  and 
the  transmission  properties  of  the  ionized  layers  which  surround 
the  earth.  The  observed  resonant  frequencies  are  presumably  characteristic 
of  the  transmission  media  and,  in  the  absence  of  a  detailed  knowledge  of 
the  sources,  provide  fruitful  grounds  for  theoretical  study. 

The  transmission  of  hydromagnetic  waves  through  regions  of 
the  ionosphere  and  exosphere  has  been  studied  by  various  workers. 

(In  what  follows  the  term  ionosphere  is  used  ;.n  connection  with  the 
altitude  region  between  about  80  km  and  500  km.)  Francis  and  Karplus 
(1960)  computed  ionospheric  magnetic  transmission  coefficients  by 
numerically  integrating  the  relevant  differential  equations,  but 
their  results  did  not  exhibit  resonances.  As  discussed  elsewhere 
(Greifinger  and  Greifinger,  1965) ,  this  was  because  they  did  not 
consider  sufficiently  low  angular  frequencies.  Jacobs  and  Watanabe 
(1962)  considered  the  transmission  of  hydromagnetic  waves  through 
the  lower  exosphere  (the  height  region  between  about  500  km  and 
2000  km),  which  they  assumed  collisionless.  By  dividing  the  lower 
exosphere  into  six  uniform  layers,  and  by  assuming  it  to  be  grounded, 
they  computed  magnetic  transmission  coefficients  which  exhibit  pronounced 
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resonances  in  the  Pc  1  and  2  range  (we  use  the  recently  suggested 
notation  of  Jaccbs  ct  al. .  (1964)).  The  lower  exosphere  has  also 
been  analyzed  by  Elliott  and  Hodder  (1965)  who  used  an  analysis 
technique  based  on  the  string  analogy.  Prince  and  Bostick  (1964), 
using  a  realistic,  but  grounded,  ionosphere  and  exosphere  have 
computed  transmission  coefficients  which  exhibit  resonances  falling 
mainly  in  the  Pc  1  and  2  categories.  Their  method  was  to  divide 
the  ionized  regions  into  a  large  number  of  uniform  layers  and  to 
proceed  numerically. 

The  purpose  of  the  present  paper  is  to  consider  theoretically 
the  transmission  of  micropulsations  through  the  ionosphere  and  lower 
exosphere.  Here,  however,  in  addition  to  the  approximations  usually 
used  to  describe  hydromagnetic  wave  propagation  in  the  ionosphere 
and  exosphere,  we  consider  the  case  of  vertical  polar  propagation 
and  represent  the  ionosphere  by  a  constant  Alfven  speed  and  locally 
exponential  ion  collision  frequency,  and  the  lower  exosphere  by 
a  collisionless  region  with  an  Alfven  speed  which  increases  exponentially 
with  height.  The  region  above  the  lower  exosphere  is  treated  as 
having  a  uniform  Alfven  speed.  By  using  this  model,  the  validity 
of  which  is  discussed  in  Section  II,  it  is  possible  to  cast  the 
ionsopheric  and  exospheric  wave  equations  into  forms  which  are 
solvable  in  terms  of  known  functions.  In  a  previous  paper  (Grei finger 
and  Grei finger,  1965)  this  was  dom  for  the  ionosphere  only,  and 
analytic  expressions  for  the  ionospheric  magnetic  transmission  and 
reflection  coefficients  were  obtained  and  numerical  results  were 
presented.  A  knowledge  of  the  contents  of  that  paper,  which  we 
will  henceforth  refer  to  at  I,  will  be  assumed  throughout  the  present 
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treatment.  Here  we  join  the  exospheric  and  ionospheric  solutions  and 
obtair  analytic  expressions  for  the  composite  transmission  and  reflection 
coefficients  which  characterize  the  entire  region  below  2000  km  or 
so.  This  procedure,  which  can  be  applied  for  oscillation  periods  in 
the  range  of  about  100  sec  to  1/10  sec,  is  described  in  Section  Ill. 

The  resulting  expressions,  although  quite  complicated,  are  less 
unwieldy  than  a  pure  numerical  treatment.  Further,  the  analytic 
forms  which  we  obtain  afford  clear  physical  interpretations  as  to 
the  origin  of  the  experimentally  measured  transmission  resonances. 

By  using  straight  forward  limiting  processes  one  can  obtain  relatively 
simple  expressions  which  describe  the  above  mentioned  cases  of 
Jacobs  and  Watanabe  (grounded  exosphere)  and  of  I  (uniform  Alfven 
speed).  These  and  other  limits  are  discussed  in  Section  IV.  In 
Section  V  we  present  numerical  results  for  several  models.  There 
it  is  shown  that,  in  addition  to  the  daytime  ionospheric  Pc  4 
resonance  (previously  called  Pc)  found  in  I  and  attributed  to  the 
filtering  effect  of  the  ionosphere,  a  series  of  higher  frequency 
exospheric  resonances  appears.  These  are  essentially  the  Pc  1  and  2 
resonances  of  Jacobs  and  Watanabe,  albeit  greatly  modified  by  the 
more  realistic  lower  boundary  conditions  used  here.  It  is  shown 
further  in  Section  V  that  the  grounded  ionosphere  model  used,  for 
example,  by  Prince  and  Bostick  (1964),  could  result  in  a  severe 
misplacing  of  the  daytime  Pc  4  resonance.  Although  the  present 
treatment  is  strictly  applicable  only  to  vertical  polar  propagation 
(it  appears  that  this  condition  could  be  relaxed)  the  computed 
transmission  coefficients  agree  satisfactorily  with  magnetic  power 
spectra  measured  at  mid  and  lew  latitudes.  This  comparison  of  theory 
with  experiment  is  made  in  Section  V. 
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II.  THE  PHYSICAL  MODEL 


As  was  shown  in  I  and  as  will  be  seen  in  the  following  Section, 
L..e  parameters  which  govern  the  vertical  polar  propagation  of 
hydromagnet ic  waves  in  the  ionosphere  are  the  ion  radial  cyclotron 
frequency,  the  ion-neut:  al  collision  frequency,  and  the 
Alfven  speed, 


c  0), 


v/4rTNieMi 


(1) 


where  c  is  the  vacuum  speed  of  light,  is  the  ionic  plasma  frequency, 

|Bq|  is  the  terrestrial  magnetic  field  strength,  is  the  volume 

density  of  ion-pairs,  e  is  L'ne  electronic  charge  (we  assume  single 

ionization),  and  is  the  average  ionic  mass  (typically  25  or  so 

proton  masses).  In  the  lower  exosphere  the  propagation  is  determined 

solely  by  the  behavior  of  V^,  provided  that  we  restrict  ourselves 

to  periods  sufficiently  short  that  collisional  effects  are  small  at 

heights  in  excess  of  500  km  or  so.  This  latter  requirement  restricts 

our  treatment  to  periods  shorter  than  about  100  sec.  It  is  quite 

reasonable  to  assume  <d^  to  be  constant  in  the  ionosphere.  Thus,  we 

wish  to  assign  ionospheric  height  dependences  to  and  (or, 

"1/2 

equivalently,  )  and  an  exospheric  height  dependence  to 

which,  in  addition  to  describing  adequately  the  physical  conditions 
of  interest,  render  the  relevant  wave  equations  solvable.  Fortunately, 
currently  used  ionospheric  and  exospheric  models  can  be  so  represented. 

In  what  follows  we  will  take  the  height  profiles  of  Alfven 
speed  used  by  Jacobs  and  Watanabe  (1962)  or  Prince  and  Bostick  (1964) 


as  being  representative.  There,  is  seen  to  be  relatively  constant 
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between  the  bottom  of  the  ionosphere  (some  80  lan  above  ground  level) 
and  (depending  upon  the  local  time  and  sunspot  conditions)  an  altitude 
of  some  4:0-500  km.  We  take  this  latter  height  as  being  the  "top" 
of  the  ionosphere.  In  the  lower  exosphere,  which  we  take  as  the 
region  between  the  ionosphere  and  seme  1500-2000  lan  (again  depending 
on  the  time  and  sunspot  activity),  exhibits  a  dependence  which  is 
remarkably  well  fitted  by  a  height-increasing  exponential.  Above 
the  lower  exosphere  V  again  becomes  quite  uniform  and  will  be 

A 

taken  as  constant.  With  regard  to  the  ion  collision  frequency, 
it  was  shown  in  I  that  by  taking  the  origin  a*-  that  height  where 

f |>  and  choosing  the  coordinate  system  so  increasing  Z  corresponds 
to  decreasing  altitude,  we  may  use 


where  p  is  a  local  scale  height,  of  the  order  of  10  km.  In  order 
for  this  representation  to  be  valid,  it  is  only  necessary  for  to 
have  the  form  (2)  within  a  few  tens  of  kilometers  on  either  side 
of  the  origin.  This  indeed  corresponds  to  the  actual  physical 
situation. 

For  *  175  rad/sec,  which  is  the  value  wh<  ch  will  be  used 
here,  the  representation  (2)  places  the  origin  some  130  km  above  the 
ground.  Thus,  our  propagation  model,  which  is  shown  schematically  in 
Fig.  1,  is  described  by  the  following  formulae: 

Z  ■  130  lan  perfectly  reflecting  ground 


130  2  Z  2:  50  km 


Earth-ionosphere  cavity 


ionosphere 


(3) 


jO  km  2i  Z  2?  Z, 


V  -  V 
VA  0 


Z2  2  z3 


V  3  V  6 
VA  V2 


(Z-Z 2) 

X 


lower  exosphere  (4) 


z  £  z. 


v  a  V  =  V  G 
VA  V3  2 


exosphere  (5) 


and  by  Eq.  (2).  The  behavior  of  at  the  level  Z^  depends  upon  the 
local  time.  The  data  indicates  that  under  daytime  conditions  Vq  = 
and  is  continuous,  but  has  a  discontinuous  derivative.  At  night, 
however,  the  depletion  of  ionization  in  the  ionosphere  causes  the 
condition  Vq  »  V£  to  obtain,  and  the  ionosphere-exosphere  boundary 
is  best  characterized  by  a  rather  severe  discontinuity  in  Alfven  speed. 
Apparently,  at  the  level  Z^i  it  is  reasonable  to  take  continuous 
with  a  discontinuous  derivative  at  all  times,  as  is  assumed  in  Eq.  (5). 
The  numerical  values  which  we  shall  use  for  these  parameters,  for 
various  local  times  and  sunspot  conditions,  will  be  tabulated  in 
the  Section  containing  our  numerical  results.  We  note  that  the 
propagation  model  shown  in  Fig.  1  is  an  extension  of  that  used  in  I, 
where  here  we  have  included  regions  above  the  ionosphere. 

In  addition  to  the  above  analytic  model  for  and  we  shall 
employ  the  same  simplifying  assumptions  as  did  Francis  and  Karplus 
(1960)  a-i  as  was  done  in  I,  viz.: 

i)  Hydrostatic  restoring  forces  and  gravitational  forces 
are  small  compared  to  electrostatic  restoring  forces  and 
may  be  neglected. 

ii)  Collisions  between  electrons  and  ions  and  collisions  of 
neutrals  with  charged  particles  may  be  neglected.  The 
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latcer  implies  that  the  neutrals  take  no  part  in  the 
mass  motion. 

iii)  The  ionosphere  is  still  and  flat, 
iv)  The  wave  normal  is  along  the  ionospheric  density  gradient. 

These  assumptions,  coupled  with  the  propagation  model  shown  in  Fig.  1, 
will  permit  an  analytic  representation  of  the  micropulsation  propagation. 
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III.  THE  COMPOSITE  TRANSMISSION  AND  REFLECTION  COEFFICIENTS 


In  paper  I  the  propagation  of  hydromagnet ic  waves  in  and 
below  the  ionosphere  was  considered,  and  ionospheric  magnetic 
transmission  and  reflection  coefficients  were  derived  and  evaluated. 

Here  we  shall  extend  this  previous  work  to  include  the  filtering 
action  of  the  lower  exosphere.  Although  a  knowledge  of  I  is  assumed 
we  shall,  for  completeness,  devote  the  first  part  of  this  Section  to 
a  very  brief  recapitulation  of  the  ionospheric  portion  of  the 
problem.  We  shall  then  go  on  to  solve  the  wave  equation  in  the 
exosphere  and,  by  joining  the  ionospheric  and  exospheric  solutions, 
shall  obtain  composite  transmission  and  reflection  coefficients. 

Since  we  are  considering  vertical  polar  propagation,  it  is 
convenient  to  introduce  the  variables 

E  -  E  -  iE  ,  B  -  B  -  iB  (6) 

x  y  x  y 

which,  for  a  time  dependence  e  *"  with  a)  >  0,  describe  the  electric 
and  magnetic  fields  of  the  left-hand  circularly  polarized  wave  (the 
L-mode).  As  discussed  in  I,  the  final  results  can  be  made  to  apply 
to  the  R-mode  simply  by  using  negative  Further,  for  hydromagnetic 
waves  in  the  ionosphere, 

(i)  od  «  I  «  |o>e| 

(li)  i>e  «  |<Del 

(ill)  Me»e  « 
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where  the  subscript  e  denotes  electron  quantities.  By  using 
dimensionless  variables  (denoted  by  a  tilde'  with  |co^  |  being  the 
unit  of  frequency  and  V  /  |o^  f  the  unit  of  length,  it  was  shown  in 
I  that  the  electric  field  satisfies  the  following  wave  equation  in 
the  ionosphere: 


where 


+  f(£,  Z/p") 


E  *  0 


(7) 


f<£,  z/<0 


l-le-Z/P 


(8) 


The  form  of  (8)  indicates  that  the  ionosphere  may  be  thought  of  as 
being  composed  of  three  regions  which  are  defined  as  follows  (these 
ionospheric  regions  should  not  be  confused  with  the  layers  shown  in 
Fig.  1): 


~  Z/n'  ~ 

i)  The  Hall  region,  where  u.  »  1  (e  w  »  1);  f  %  od  an  l 

1  +i£*  I 

is  a  linear  combination  of  the  solutions  e~ 

***  ^ 

ii)  The  Alfven  region,  where  v./(5  «  1  (e^^  «  59 ;  f  s»  7Z? 
and  E  is  a  combination  of  the  solutions  e~  .  Mote 
that  for  very  long  periods,  say  cn  £  .05  rad/sec,  there 
is  no  Alfven  region  below  the  lower  exosphere.  This 
will  place  a  limitation  upon  the  forthcoming  analysis, 

iii)  The  transition  region,  which  falls  between  the  above  two 
regions,  and  where  no  simplification  of  Eqs.  (7)  and  (8) 
can  be  made. 

Equation  (7)  may  be  transformed  to  a  form  of  the  hypergeometric 

■z/y 


equati.»T'  with  u  ■  1  -  it 


as  an  independent  variable.  Thus,  the 
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transition  region  elution  for  £  involves  a  superposition  of  hyper¬ 
geometric  functions,  which  goes  to  the  proper  limiting  form  in  the 
Hall  and  Alfven  regions.  This,  in  fact,  provides  the  connection 
between  the  fields  in  the  Alfven  region  and  those  at  the  bottom  of 
the  ionosphere.  By  carrying  out  the  mathematics  of  this  connection, 
and  by  applying  the  boundary  condition 


at  the  bottom  of  the  Hall  region,  it  was  shown  in  I  that  the  fields 
in  the  Alfven  region  may  be  written 


^  E(Z.  ,uj)  . * — rx 

E(Z)  -  - J -  (e^Z+  re~UiiZ) 

E 


B(Z)  = 


BCZ^S) 


(ID 


where  t  and  t  are  the  electric  and  magnetic  transmission  coefficients 
£  B 

of  the  ionosphere,  i.e.,  the  ratio  of  the  total  fields  at  Z  =  Z^  to 

the  fields  of  the  downcaning  part  of  the  signal  in  the  Alfven  region. 

Clearly,  r  is  the  ionospheric  reflection  coefficient.  The  forms 

of  t  ,  t  ,  and  r  were  found  in  I  to  be 
£  B 


-li¬ 


en) 


(14) 

and  r  denotes  the  Gamma-function.  Also,  so  far  as  the  functional 
form  of  g  is  concerned,  S  and  S  are  treated  as  independent. 

We  now  proceed  to  the  main  point  of  the  present  paper,  namely 
'to  find  solutions  to  the  wave  equation  in  the  exosphere  and  to  join 
these  solutions  with  those  just  obtained  for  the  ionosphere.  For 
the  purpose  of  this  joining,  we  shall  consider  oscillation  periods 
which  are  sufficiently  short  (—  £  100  sec)  that  E  and  B  have  the 

<D 

limiting  forms  (10)  and  (11)  near  the  top  of  the  ionosphere.  For 
much  longer  periods  this  joining  would  be  quite  complicated,  since 
the  transition  region  would  then  extend  up  to  or  beyond  the  ionosphere - 
exosphere  boundary,  and  the  fields  just  beneath  this  boundary  would 
have  to  be  expressed  in  terms  of  hypergeometric  functions. 

According  to  our  model  shown  in  Fig.  1,  the  medium  is  uniform 
for  Z  <  Z^  and  the  hydromagnetic  wave  equation,  subject  to  the 
above  mentioned  restrictions  on  frequency,  is 
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We  may  write 


(16) 


(17) 


for  Z  s  Z_,  where  T  ,T  and  R  are  the  composite  transmission  anu  reflection 
3  BE 

coefficients  for  the  entire  region  Z^  >  Z  ^  Z^,  and  are  to  be  found. 

The  application  of  the  induction  equation  at  Z^  and  at  Z  <  Z^  shows 


that 


(18) 


In  order  to  connect  the  solutions  (16'  and  (17)  to  the 
solutions  (10)  and  (11),  we  must  find  the  form  of  the  wave  fields 
in  the  lover  exosphere.  Again  referring  to  our  model,  the  wave 
equation  in  the  rev  jn  Z^  Z  ^  Z^  is 


(19) 


(20) 


as  an  independent  variable  becomes 
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(  '•-2  d  ~  d  —2  — 2  ~2  \ 

y  T^+y~  +  "'  x  y  e 

\  dy  dy  j 


Equation  (21)  is  of  course  Bessel's  equation  of  zero  order,  and  the 
electric  field  in  the  lower  exosphere  may  be  written 


E  =  A(oj)(J  (S  \  y)  +  D(t5)  N  (u>  X  y))  Z0  >  Z  ^  Z, 


where  J  and  N  are  the  Bessel  and  Neumann  functions,  and  A  and  D 
o  o 

are  to  be  determined  from  an  application  of  the  boundary  conditions . 

^  .Srf 

Now,  by  requiring  that  E  and  dE/d7.  be  continuous  at  Z  »  Z^  and 
Z  *  Z  ,  and  by  using  the  forms  (10),  (16)  and  (22),  four  readily 
solvable  equations  for  A,D,R,  and  T  result.  By  solving  these,  and 

£i 

by  using  Eq.  (18),  we  find  the  quantities  of  interest,  T  and  R,  to  be 

B 

-  z,) 

V 

~  3 

4iV  e  J 

T  =  3.  .  _  ,  .  —  ■  — 

B  ~  r 

n  oj  X 


81  Jq((o  X  y3)  +  iJ^(to  X  y3)  J  -  or  I  Nq(o)  X  y3>  +  iN1(to  X  y3)  | 


where 


O'  ■  JQ(a>  X  y2)  -  iy2J1(a)  X  y2)  -  re 


■2ia>Z, 


|^Jo(cb  X  y2)+i^J1(co  K$]  (24) 
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2ia>Z, 


3  =  Nq((o  \  y2)  -  iy^CS  x  y2>  -  re  2  NQ(m  \  y?)  +  iy^N^co  X  (25) 


2i  f-  Z„ 

\  V  3/J  (m  \  yj  +  GN  (Z  X  yj  -  i[J.(ro  x  yJ  +  GN  (S  r  yJ] 

R  =  e  3  -2 - 2 - 2 - - - - - 2 - - - 2  —  (26) 

Jo(co  \  y3)  J-  GN  (oj  x’  y3)  +  i[J1(o>  X  y3)  +  GN^cd  X  y3)] 


where 


ro>Z0  -ituZ  ^  _ 

e  [JQ(S  X  y2)-ly2J1(®  ~  72)>re  ^-J0^  ~  y2^+iy2Jl^  *  3^ 

G  *  -  - •  (27) 

„  “UoZ,} 

e  [Nq(o)  X  y2>‘-iy2N1(cD  X  y2>  >re  [N0(®  X  y^+iy^^cD  X  ^)] 


In  these  expressions,  t  and  r  are  the  ionospheric  transmission  and 

U 

reflection  coefficients  which  were  defined  earlier,  and  jr  means 
'y,  as  defined  by  Eq.  (20),  evaluated  at  Z  =  Z2>  etc.  Finally,  the 
Wronskian  relation 


J,  (cd  X  y)  N  (a)  x  y)  -  J  (cd  x  y)  N,(a>  x  y)  s  - 

1  O  O  JL  ^  ~ 

n  cd  X  y 


was  used  in  deriving  Eqs.  ( 23) - ( 27 ) - 
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IV.  SOME  LIMITING  CASES 

The  expressions  (23)  and  (26)  for  T  and  R  have  been  numerically 

evaluated  for  several  propagation  models,  and  the  results  will  be 

given  in  the  following  Section.  Prior  to  this,  however,  we  shall 

consider  some  of  the  limiting  forms  of  T  and  R.  These  limiting 

B 

cases,  which  are  much  simpler  than  the  general  expressions  shown 
above,  will  prove  very  useful  in  interpreting  the  forthcoming 
numerical  results.  Further,  by  considering  limiting  situations,  it 
will  be  possible  to  draw  conclusions  regarding  the  origin  of  many 
of  the  experimentally  measured  transmission  resonances. 

Jacobs  and  Watanabe  (1962)  considered  the  transmission  of 
hydromagnet ic  waves  through  a  grounded  lower  exosphere,  and  obtained 
pronounced  transmission  resonances.  Their  model,  is,  in  fact,  a 
limiting  case  of  that  shown  in  Fig.  1.  We  first  replace  the  ionosphere 
of  our  model  by  a  perfectly  conducting  ground,  which  is  accomplished 
by  setting  p  »  H  *  Zj=  0.  This  leads  to  the  values 

r  *  -1;  *  2  (29) 

for  the  ionospheric  coefficients,  while  a  and  &  become  simply 

a  M  2JQ((o\y2);  P3  2Nq(S  X  y2)  .  (30) 

At  the  same  time,  we  must  also  preserve  the  distance  between  Z^  and  Z y 
Thus,  since  we  are  lowering  the  level  of  the  bottom  of  the  exosphere 
by  an  amount  j | >  we  must  also  replace  Z^  by  -  ?  +  Z.^  (remembering 

that  Z^  <  0).  Making  these  changes  in  Eq.  (23),  we  see  that  the 
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magnetic  transmiss '.on  coefficient  for  the  grounded  exosphere  is 
given  by 


UtV}e 


Grounded 

Exosphere 


TT  03  X 


[No(u>  X  y2)  [JoG  X  y^)+iJ1(tP  X  yy } 


JoGx72)  iN0(Sr73)+iN1®r73)}]'1  (31) 


Wi.ere  y^  and  y^  are  now  referred  to  the  new  co-ordinate  system. 

As  will  be  seen  presently,  the  magnitude  of  T  as  given  by  Eq.  (31) 
agrees  precisely  with  the  numerical  results  of  Jacobs  and  Watanabe. 

The  transmission  resonances,  which  occur  at  the  minima  of  the  denominator, 
are  associated  solely  with  the  filtering  ac  jn  of  the  lower  exosphere. 

We  can  deduce  very  simple  relations  for  the  location  and  magnitudes  of 
these  resonances  by  further  simplifying  Eq.  (31).  The  numerical 
values  for  the  parameters  used  in  the  model  will  be  given  in  the  next 
section,  but  it  turns  out  that,  generally  speaking, 


y 


3 


(32) 


and  ,  as  will  be  seen  a  posceriori .  the  first  few  exospheric  resonances 
are  so  located  that 


a)  \  y^  «  !• 


(33) 
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Thus , 


Jo(cd  x  y3)  +  UjCa)  \  >’3)  ~  1 


(34) 


|No(ST  y3)  +  X  y3)  I  ~  |“  In  (  2 


-)  +  - - - 1  »  1  (  35) 

^  .i^ 


1.8  £1)  )  y3  TT  OD  X  y3 


and  since  N  (oo  X  y_)  is  of  order  unity  or  less,  the  first  few  grounded 
o  2 

exospheric  resonances  are  located  such  that 


Jo(tt>  X  y2)  0 


(36) 


l.  e.  , 


O) 


qn 
X  y. 


(37) 


*  th 

where  q^  is  the  n  root  of  Jq  (q^  *  2.4,  q^  “  5.5,  q3  =  8.6  •••). 
The  combination  of  Eqs.  (32)  and  (37)  justifies  the  inequality  (33). 
The  roots  of  Jq  may  be  quite  accurately  located  by  using  the 


asymptotic  form  J  (x)  ~  cos(x-tt/4).  Now  reverting  to 

o  ttx 

conventional  units,  we  note  that,  in  view  of  Eqs.  (20)  and  (32),  the 
phase  change  Acp  between  and  is 


5cp  =  f  dZ  =  Xa> 

7  A  4*0  *  9 J 


Z3  ’A 


Xflo 

V„ 


(38) 


which  is  just  the  argument  of  Jq  in  Eq.  (36)  Thus  the  location  of  the 
resonant  frequencies  can  be  written  approximately 


rr 

cos  ( Acp  -  7)  «  0 
n  4 


(39) 


so  that,  in  the  sense  of  a  phase  integral,  the  grounded  lower  exosphere 
deviates  from  a  quarter-wave  plate  behavior  only  by  the  phase 
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factor  tt/4.  The  magnitudes  of  these  resonant  peaks  are  also  quite 
easily  found.  By  using  (34),  (36)  and  (j7)  in  conjunction  with  (31) 
it  follows  that 


Grounded 

Exosphere 


4  V  y 
3y2 

rrq  N  (q  ) 
^n  o  nn 


(40) 


As  will  be  seen  presently,  when  the  exosphere  is  terminated  upon  the 
Alfven  region  of  the  ionosphere,  rather  than  upon  a  perfect  earth, 
the  transmission  peaks  are  shifted,  broadened  and  depressed. 

Another  limit  which  provides  useful  insight  is  that  of  the 
uniform  exosphere.  This  case,  which  is  mathematically  represented 
by  taking  \  *  oa}  and  hence  y^  **  y^  and  =  removes  the  filtering  action 
of  the  lower  exosphere  and  thus  eliminates  the  exospheric  resonances. 

In  this  case,  by  using  the  Wronskian  relation  (28),  Eqs.  (23)  and 
(26)  become 


T 


2S 

l+y2+re  (y2"l) 


(41) 


and 


21  Z3  _  ~2UoZ  -2i^Z 

3  y,(l  +  re  Z)  -  (1  -  re  ) 

- 2 - .  (42) 

_  -2&Z  -  2  itnzl 

y2(l  +  re  )  +  (1  -  re  ) 


Since  the  exospheric  resonances  have  been  removed  in  this 
limit,  any  resonances  which  remain  must  obviously  be  associated  with 
the  filtering  action  of  the  ionosphere.  For  daytime  conditions, 
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>2  °  1  and  ■  1,  so  that  Eqs.  (41)  and  (42)  become 


y 


<■“  -j 

1<d(Z  -Z  ) 

J  z  t 

B 


(43) 


(44) 


These  are,  of  course,  just  the  ionospheric  transmission  and  reflection 
coefficients  referred  to  th*  height  Z^,  and  the  resonances  are 
consequently  exactly  those  obtained  in  I.  On  the  other  hand,  in  the 
nighttime  a  discontinuity  in  Alfven  speed,  and  hence  an  impedance 
mismatch,  arises  at  the  height  Z^.  The  result  of  this  mismatch,  which 
is  rather  severe,  is  the  virtual  elimination  of  the  lowest  (and,  for 
practical  purposes,  only)  ionospheric  resonance. 

The  numerical  results  for  these  limits,  together  with  those  for 
the  full  expression,  will  be  presented  in  the  next  section. 
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V.  NUMERICAL  RESULTS  AND  DISCUSSION 


In  this  Section  we  presert  numerical  results  for  the  transmission 
and  reflection  coefficients  for  several  model  ionospheres  and  lower 
exospheres.  These  results,  although  obtained  for  polar  propagation, 
will  be  compared  with,  and  will  be  seen  to  agree  well  with, 
experimental  data  obtained  at  mid  and  low  latitudes.  The  relevant 
parameters,  which  were  chosen  essentially  to  correspond  to  the 
Alfven-speed  height-profiles  used  by  Jacobs  and  Watanabe  (1962)  or 
Prince  and  Bostick  (1964)  are  shown  in  Table  1.  In  addition,  we 
shall  use  |<jd^|  =  175  rad/sec  and  an  ion-collision  frequency  scale 
height,  p,  of  10  km.  These  two  values  were  shown  in  I  to 
describe  conditions  in  the  ionosphere  satisfactorily. 


Local  time  and  sunspot 
conditions 

V  (km/sec) 
o 

^^(lcm/  sec) 

\(km) 

Z2(km) 

Z3(km) 

Day-Sunspot  Max. 

200 

200 

445 

-220 

-1870 

Night-Sunspot  Max. 

5000 

450 

455 

-370 

-1950 

Day-Sunspot  Min. 

400 

400 

310 

-320 

-1320 

Night-Sunspot  Min. 

5000 

800 

280 

-270 

-1170 

Table  1 

IONOSPHERIC  AND  LOWER  EXOSPHERIC  PARAMETERS 
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The  quantity  of  interest  with  respect  to  understanding  the 

observed  spectra  of  micropulsations  is  the  composite  transmission 

coefficient.  Our  numerical  results  for  the  magnitude  of  T_  for  the 

B 

L-mode,  corresponding  to  the  four  models  given  in  Table  1,  are 
shown  in  Figs.  2-5.  The  R-mode  transmission  coefficients,  for  the 
day  and  night  sunspot  minimum  models,  are  shown  in  Figs.  6  and  7. 

We  have  also  shown  in  Figs.  2-7,  for  purposes  of  comparison,  the 
ionospheric  transmission  coefficients  | t^ j  and  | (X  =  ®) | .  It  will  be 
recalled  that  for  our  daytime  models,  these  two  quantities  become 
identical.  Also,  we  have  included  in  Fig.  2  the  results  obtained 
by  taking  the  grounded  exosphere  limit  considered  by  Jacobs  and 
Watanabe  (1962)  and  represented  analytically  by  our  Eq.  (31). 

The  L-mode  transmission  coefficient,  computed  for  daytime  sunspot 
maximum  conditions,  but  with  H  *  0  (grounded  ionosphere)  is  shown 
in  Fig.  8.  The  composite  reflection  coefficients  are  of  limited 
use  in  interpreting  experimental  observations,  and  exhibit  a  generally 
complicated  structure  of  sharp,  narrow  minima  which  occur  near  the 
transmission  resonances.  Sample  curves  of  R  and  r,  for  daytime 
sunspot  minimum  conditions,  are  shown  in  Fig.  9. 

We  start  by  discussing  the  daytime  cases  and  note,  by  comparing 
Figs.  2  and  3  with  Fig.  6,  that  the  daytime  resonant  structure  of 
the  transmitted  signal  is  governed  primarily  by  the  L-mcde.  This 
would  imply,  at  least  for  high  latitudes,  that  the  polarization  of 
the  daytime  signals  near  the  resonant  frequencies  should  be  nearly 
circular.  The  results  shown  in  Figs.  2  and  3  indicate  that  the 
main  conclusion  drawn  in  I  is  unaffected  by  the  inclusion  of  the 
lower  exosphere.  In  particular  the  daytime  Pc  4  resonance,  which  has 
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a  period  of  many  tens  of  seconds  (a  &  10  "*)  and  which  is  clearly 
due  to  the  filtering  action  of  the  ionosphere,  remains  intact.  Still 
considering  the  two  daytime  models,  we  see  by  comparing  It  I  with 
|Tgj  in  Figs.  2  and  3  that  the  inclusion  of  the  lower  exosphere 
introduces  a  series  of  higher  frequency  resonances  falling  in  the 
lower  Pc-3  to  Pc  1  range.  These  peaks,  which  are  clearly  exospheric 
in  origin,  are  a  modified  version  of  those  predicted  by  Jacobs  and 
Watanabe  (1962).  To  illustrate  this  point  we  have  included  a  plot 
of  |Tg|  for  a  grounded  exosphere  (see  Eq.  (31))  in  Fig.  2.  This 
curve,  which  agrees  closely  with  the  corresponding  numerical  results 
of  Jacobs  and  Watanabe,  exhibits  resonances  which  are  larger,  more 
narrow,  and  shifted  in  frequency  from  those  evident  in  |T^|  as 
obtained  for  the  more  accurate  model. 

One  of  the  most  striking  features  of  the  daytime  models  is 
that  the  first  exospheric  resonance  is  comparable  in  magnitude  with 
that  due  to  the  ionosphere.  This  pronounced  double-peak  in  the 
daytime  spectrum  is  in  good  qualitative  agreement  with  the 
two  prominent  daytime  bands  measured  by  several  workers.  Quanti¬ 
tatively,  our  computed  values  for  the  location  and  spacing  of 
these  two  peaks  agree  better  with  the  measurements  of  Ness ,  et  al.,  (1962) 
and  Smith,  et  al. ,  (1961)  than  with  those  of  Maple  (1959)  and  Davidson 
(1964).  Here  we  see  that  the  two  peaks  arise  in  different  ionized  layers. 
There  is  also  a  series  of  higher  frequency  exospheric  Pc  1  and  2 
daytime  resonances.  The  magnitude  of  these  falls  off  rapidly  with 
increasing  frequency,  due  to  the  opacity  of  the  daytime  ionosphere. 

It  is  thus  questionable  if  any  but  perhaos  the  first  few  cf  these 
could  be  easily  measured. 
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Prince  and  Bostick  (1964)  have  also  computed  magnetic  transmission 
coefficients,  which,  over  much  of  the  spectrum,  appear  similar  to 
those  shown  here.  A  major  difference,  however,  is  the  absence  of 
the  experimentally  measured  Pc  4  daytime  resonance  in  their  work. 

Although  this  could  possibly  be  a  consequence  of  the  fact  that  their 
vertical  equatorial  propagation  model  was  different  than  the  polar 
situation  treated  here,  we  feel  that  another  likely  explanation  is 
their  use  of  a  grounded  ionosphere  (H  *  0).  We  have  shown,  in 
Fig.  8,  |Tg|  computed  for  the  same  daytime  maximum  conditions  as 
in  Fig.  2,  except  that  H  has  been  taken  equal  to  zero.  A  comparison 
of  these  two  figures  shows  that  the  assumption  of  a  grounded  ionosphere 
causes  the  first  daytime  resonance  to  be  depressed  and  shifted  to 
a  considerably  higher  frequency.  The  cause  of  this  behavior  was 
discussed  in  I.  The  fact  that  the  transmission  coefficients  computed 
by  Prince  and  Bostick  fall  off  with  increasing  frequency  much  faster 
than  those  shown  here  is  due  to  the  differences  between  polar  and 
equatorial  propagation. 

We  now  discuss  the  nighttime  transmission  coefficients,  which 

differ  both  qualitatively  and  quantitatively  from  those  obtained 

for  daytime  conditions.  One  major  difference,  which  can  be  seen 

in  Figs.  4  and  5,  is  that  the  L-mode  ionospheric  Pc  4  resonance, 

which  was  prominent  in  the  daytime,  essentially  vanishes.  Since 

neither  t  nor  T  (\  **  <=)  (which  includes  the  effects  of  the 
B  B 

discontinuity  in  V.  at  Z  ■  Z  )  exhibiLs  any  pronounced  nighttime 

A  4 

resonances,  the  resonant  peaks  in  T  are  due  to  the  lower  exosphere. 

These  exospheric  resonances  are  much  more  pronounced  at  night  than 
in  the  daytime.  Finally,  and  of  equal  importance,  Fig.  7 
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shows  that  the  first  few  R-mode  resonances  are  nearly  identical  in 
size  and  location  to  those  for  the  L-mode.  This  is  a  consequence 
of  the  relative  transparency  of  the  nighttime  ionosphere.  Thus, 
while  we  would  expect  high  latitude  daytime  signals  to  be  nearly 
circularly  polarized  (at  least  at  the  resonant  frequencies)  nighttime 
ones  would  have  essentially  the  polarization  of  the  sources  (assuming 
a  steady  state).  Detailed  measurements  of  the  resonant- fiequency 
polarization  of  high- latitude  magnetic  signals  would  be  a  useful 
check  on  the  results  presented  here. 

The  resonant  structure  of  the  nighttime  transmission  coefficients 
agrees  at  least  qualitatively  with  many  of  the  experii  ^ntal  data.  For 
example,  we  interpret  the  first  and  largest  nighttime  transmission 
peak  as  the  single  prominent  nighttime  band  noted  by  Maple  (1959) 
although  the  resonant  frequencies  computed  here  are  perhaps  twice 
as  large  as  the  1/8  cps  which  he  measured.  Also,  the  computed 
nighttime  increase  in  the  magnitude  of  the  Fc  1  and  2  resonances 
is  in  agreement  with  the  nocturnal  appearance  of  certain  fluctuations 
in  this  frequency  range  (Feniof f ,  1960;  Tepley  and  Amundsen,  1965). 

We  must  emphasize  that  the  diurnal  variation  of  the  sources  could  be 
an  important  factor  here,  however.  Thi  possibility  of  strong 
daytime  sources  overcoming  the  generally  poorer  daytime  Pc  1  and  2 
transmission  must,  of  course,  also  be  considered. 

In  conclusion,  our  analytic,  polar  propagation  model  has 
yielded  readily  interpretable  results  which  agree  with  experiment 
in  the  following  respects: 

1)  A  strong  double  daytime  resonance  in  agreement  with 
the  measurements  of  Maple  (1959),  and  others  has  been 
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2)  A  strong  nighttime  resonance,  again  in  agreement  with 
Maple ,  has  been  obtained. 

3)  The  computed  nocLurnal  prominence  of  the  exospheric 

Pc  1  and  2  resonances  agrees  with  the  measured  diurnal 
behavior  of  certain  classes  of  micropulsations. 

The  basic  structure  of  the  micropulsation  spectra  is  thus  contained 
in  the  expressions  given  here.  A  more  detailed  quantitative  agreement 
with  experiment  could  presumably  be  attained  by  adjusting  the  values 
of  the  parameters  used.  Further,  the  ability  to  take  simple  limiting 
cases  of  our  analytic  expressions  has  provided  physical  insight  with 
regard  to  the  transmission  and  reflection  processes.  There  are,  of 
course,  limitations  in  our  treatment.  The  most  severe  is  perhaps  the 
restriction  to  polar  propagation.  This  prevents  an  analysis  of 
the  latitude  dependence  of  the  transmission  resonances  and  polarization 
of  ground-level  signals.  The  extension  to  non-polar  situations 
does  seem  feasible,  however,  and  is  presently  being  considered. 
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Fig.  8 
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